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A3STRACT  Three  test  series  of  single-edge  notched  beams  in  three-point  bending  were 
conducted  to  evaluate  the  fracture  energy  of  concrete.  The  fracture  energy  was  determined 
from  the  area  under  the  complete  load  versus  load-point  deflection  diagram.  The  nonlinear 
Fictitious  Crack  Model  was  implemented  in  a  finite  element  analysis,  showing  good  agreement 
with  the  experimental  data. 

By  varying  the  notch  depth  and  the  beam  depth  it  was  shown  that  the  fracture  energy, 
traditionally  presented  as  a  material  property,  depends  upon  the  specimen  size  and 
configuration.  This  is  attributed  to  the  energy  dissipation  in  the  process  zone  which  is  not 
accounted  for  in  the  analytical  model. 
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’ By  varying  the  notch  depth  and  the  beam  depth  it  was  shown  that  the 
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Three  test  series  of  single-edge  notched  beams  in  three-point 
bending  were  conducted  to  evaluate  the  fracture  energy  of  concrete. 
The  fracture  energy  was  determined  from  the  area  under  the  complete 
load  versus  load-point  deflection  diagram.  The  nonlinear  Fictitious 
Crack  Model  was  implemented  in  a  finite  element  analysis,  showing 
good  agreement  with  the  experimental  data. 

By  varying  the  notch  depth  and  the  beam  depth  it  was  shown 
that  the  fracture  energy,  traditionally  presented  as  a  material 
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INTRODUCTION 


Classical  linear  elastic  fracture  mechanics  (LEFM),  which  has  been 
successfully  applied  to  metallic  and  brittle  materials,  is  limited  when 
applied  to  concrete.  LEFM  cannot  model  the  behavior  of  small  specimens 
of  the  size  typically  used  In  laboratories.  As  a  consequence,  several 
nonlinear  models  have  been  developed  which  approximate  LEFM  for  large 
sizes.  The  two  best  known  models  are  the  fictitious  crack  model  (FCM) 
introduced  by  Hillerborg  et  al.  (Ref  1),  and  the  smeared  or  crack  band 
model  (CBM)  introduced  by  Rashid  and  developed  by  Bazant  et  al .  (Ref  2, 
3,  4,  5).  Among  others,  the  two-parameter  model  by  Jenq  and  Shah 

(Ref  6)  is  more  recent  and  is  supported  by  only  limited  data. 

In  finite  element  applications,  CBM  approach  shows  a  dependency  on 
the  element  size  used  in  the  mesh.  Results  do  not  converge  in  succes¬ 
sive  analyses  where  the  element  size  is  continuously  reduced.  This 
problem  can  be  circumvented  by  linking  the  softening  stiffness  of  the 
cracked  elements  to  the  fracture  energy  G^.  However,  determined  from 
the  load-deflection  diagram  is  suspected  to  depend  on  specimen  geometry 
and  size  (Ref  5).  The  FCM  is  also  based  on  G^,  which  is  assumed  to  be  a 
material  property. 

The  existing  approaches  for  determination  of  the  fracture  energy 
are  evaluated  in  this  report  and  a  new  test  method  proposed.  Three 
series  of  tests  and  a  finite  element  analysis  using  FCM  were  conducted. 
The  objectives  were: 

•  Calculate  the  fracture  energy  following  RILEM  Technical 
Committee  50-FMC  guidelines  (Ref  7)  and  compare  it  to  the 
fracture  energy,  G^*,  according  to  a  method  presented  in 
this  report. 


•  Conduct  a  finite  element  analysis  with  different  strain 
softening  relations  based  on  the  experimental  value  of  G^. 

•  Determine  the  effect  of  notch-to-depth  ratio  and  specimen 
size  on  G^. 


EVALUATION  OF  THE  FRACTURE  ENERGY  Gf 

For  a  beam  in  three-point  bending  (Figure  la),  the  load  typically 
varies  with  load-point  deflection  (LLPD)  as  shown  in  Figure  lb. 

The  LLPD  plot  comprises  three  stages  of  behavior.  The  deflection 
increases  linearly  with  the  load  in  the  first  stage  and  the  crack,  is 
opened  but  does  not  extend.  A  fracture  process  zone  develops  during  the 
second  stage  where  microcracks  form  and  slow  crack  growth  is  apparent. 
In  the  third  stage,  known  as  the  strain  softening  zone,  rapid  crack 
growth  is  evident.  During  strain  softening  most  of  the  damage  to  the 
specimen  is  concentrated  in  a  narrow  zone.  This  concentration  is  higher 
as  the  load  carrying  capacity  decreases  and  energy  dissipation  even¬ 
tually  occurs  through  a  single  major  crack.  Strain  softening  has  been 
considered  a  material  characteristic. 

The  area,  UQ,  under  the  LLPD  curve  (Figure  2a)  represents  the 
energy  required  to  break  the  specimen.  For  a  single-edge  notched  beam 
in  three-point  bending,  RILEM  TC  50-FMC  defined  the  fracture  energy,  G^, 


Gf  =  (UQ+  mgdo)/A 


where  A  =  ligament  area  =  B(W-ao) 

B  =  width 

W  =  specimen  depth 

aQ  =  notch  depth 

mg  =  weight  of  the  specimen 

d  =  load-point  deflection  at  fracture 
o  r 
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This  definition  relies  on  the  assumption  that  all  the  energy  re¬ 
quired  to  break  the  specimen  is  transformed  into  surface  energy  by 
extension  of  a  single  macrocrack.  However,  energy  dissipation  outside 
of  the  fracture  zone  is  included  in  the  determination  of  and  should 
not  be  overlooked.  This  energy  depends  on  the  specimen  size  and  notch 
depth.  It  is  dissipated  in  creating  and  extending  a  process  zone  by 
debonding  aggregates  and  opening  microcracks.  Most  of  this  energy  con¬ 
sumption  is  believed  to  occur  during  slow  crack  growth.  Consequently, 
using  the  whole  area  under  the  curve  leads  to  an  overestimate  of  G^. 
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From  the  preceding  observations  it  was  concluded  that  an  improved 
measure  of  the  toughness  of  concrete  would  be  obtained  with  a  three- 
point  bend  beam  specimen  and  the  following  procedure: 

1.  Set  up  the  beam  with  the  notch  on  the  top  surface.  This  will 
help  in  applying  dye  into  the  cracked  surface  for  determining  crack 
growth. 

2.  Load  the  specimen  up  to  the  point  of  instability  defined  as  the 
point  past  peak  load  where  the  load  drops  ofi  to  95  percent  of  its  maxi¬ 
mum  value,  then  remove  the  load  completely.  The  area  enclosed  in  this 
load-unload  loop  includes  the  energy  spent  on  formation  of  a  process 
zone,  slow  crack  growth,  and  the  inelastic  energy  spent  outside  of  the 
crack  zone. 

3.  Insert  dye  through  the  notch  and  allow  it  to  flow  into  the 
crack  to  highlight  the  crack  length,  a  ,  at  the  point  of  instability. 

4.  Reapply  the  load  and  obtain  the  strain  softening  zone.  A  one¬ 
time  unloading  and  reloading  will  not  significantly  affect  the  LLPD 
curve. 
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5.  Define  U  as  the  total  area  under  the  LLPD  curve  minus  the  area 
in  the  load-unload  loop  Indicated  in  the  second  step  (see  Figure  2b). 

6.  Define  G^*  as  the  energy  spent  on  developing  one  major  crack 
divided  by  the  ligament  area  existing  at  that  moment: 

Gf*  -  U/B(W-ap) 

DETAILS  OF  EXPERIMENTAL  METHODOLOGY 

Fracture  Area 

During  the  strain  softening  process  a  crack  will  actually  follow  a 
surface  which  is  not  flat  but  governed  by  the  aggregate  size  and  rela¬ 
tive  hardness  compared  to  the  mortar  matrix.  An  invariant  will  be 
obtained  using  the  ligament  area,  B(W~ao),  if  the  total  energy  spent  in 
the  formation  of  a  unit  projected  area  is  constant  (as  assumed  in  the 
RILEM  approach),  or  if  the  energy  spent  other  than  in  formation  of  the 
macrocrack  is  discarded  (as  attempted  here). 

Precracking 

Precracking  of  the  specimen  (or  fatigue  cracking)  is  not  necessary 
in  the  proposed  method.  Measurement  of  the  energy  takes  place  only 
after  a  sharp  crack  has  been  formed,  and  does  not  depend  on  the  initial 
condition,  whether  it  is  precracked,  form  notched,  or  saw  notched. 

Rate  of  Loading 

RILEM  recommends  reaching  peak  load  after  30  to  60  seconds  which 

-6  -4 

corresponds  to  a  rate  on  the  order  of  5  x  10  m/sec  (2  x  10  in. /sec). 
The  work  of  fracture  and  the  strain  energy  release  rate  Increase 
slightly  (about  15  percent)  for  cross-head  deflection  rates  from  5  x 
10  ^  to  5  x  10  ^  m/sec  (2  x  10  ^  to  2  x  10  ^  in. /sec)  (Ref  8). 


Beam  Weight 


By  supporting  half  of  the  beam  weight  at  the  ends,  a  complete  LLPD 
curve  is  obtained.  The  LLPD  curve  actually  begins  after  the  applied 
load  equals  half  of  the  specimen  weight.  In  the  LLPD  plot  the  effect  of 
beam  weight  is  easily  identified  and  discarded  from  the  total  area  under 
the  curve  (Figure  6).  This  discarded  area  (1/2  weight  x  midspan  deflec¬ 
tion  at  total  fracture)  corresponds  to  work  being  transformed  into  po¬ 
tential  energy  as  the  center  of  gravity  of  the  beam  is  forced  up. 

Displacement  and  Support  Indentation 

According  to  RILEM,  the  midspan  deflection  can  be  measured  with 
reference  to  the  loading  apparatus  as  long  as  the  inelastic  indentations 
at  the  load  points  do  not  exceed  0.01  mm  (0.0004  in.)  (Ref  7).  For 
nonstandard  specimens  with  small  span/depth  ratios  (e.g.,  S/W  =  4)  the 
inelastic  indentations  at  the  loading  points  are  not  negligible.  They 
have  to  be  considered  or  else  the  deflection  due  to  load  indentations 
has  to  be  measured. 

The  error  caused  by  inelastic  indentation  is  estimated  for  a  27.65 
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MN/m  (4,000  psi)  concrete  specimen  with  dimensions  102  by  7  by  76  by 
406  mm  (4  by  3  by  16  in.)  (depth  by  width  by  span),  with  an  initial 
notch  depth  aQ  =  25  mm  (1  in.)  and  bearing  directly  on  51-mm  (2-in.) 
diameter  rollers.  A  maximum  load  of  approx imately  3.12  kN  (700  lb) 
should  be  expected.  The  minimum  bearing  area  at  the  center  roller  is 
3.12/0.02765  =  113  mm^  (0.175  in.^)  and  the  minimum  bearing  width  is 
113/76  =  1.5  mm  (0.058  in.).  This  implies  an  indentation  of  the  flat 
surface  at  the  center  roller  only  of  0.75  x  0.75/25.5  =  0,022  mm  (about 
0.001  in.).  In  this  case  the  indentation  represents  about  25  percent  of 
the  midspan  deflection  at  peak  load. 

The  clip  gage  described  by  ASTM  E399  seemed  most  appropriate  for 
accurate  displacement  measurements.  Clip  gages  were  manufactured  out  of 
high  strength  aluminum  (7075-T6)  which  was  more  readily  available  and 
easier  to  machine  than  a  titanium  alloy  as  recommended  by  ASTM.  High 
strength  aluminum  presents  a  ratio  of  yield  strength  to  modulus  of 


elasticity  as  high  as  0.0069.  High  strength  aluminum  ensures  a  large 
range  of  measurement  without  permanent  deformation  of  the  gage.  Two 
clip  gages  were  employed,  one  on  each  side  of  the  specimen,  to  mitigate 
errors  due  to  asymmetry. 

Point  of  Instability 

The  point  of  instability  was  chosen  as  the  point  after  maximum  load 
where  the  load  decreases  to  95  percent  of  its  peak  value  as  recommended 
by  Swartz  and  Yap  (Ref  9).  A  small  variation  of  load  near  peak  value  is 
accompanied  by  a  small  displacement  on  the  LLPD  curve;  however,  this 
small  amount  of  external  work  causes  a  significant  crack  advance.  This 
is  apparent  on  typical  load  versus  crack  mouth  opening  displacement 
(LCMOD)  plots  (Ref  9)  where  the  CMOD  increases  significantly  for  almost 
constant  maximum  load.  This  instability  is  attributed  to  a  redistribu¬ 
tion  of  the  elastic  energy  to  surface  energy  inside  the  specimen.  Thus, 
it  is  necessary  to  measure  crack  length  past  peak  load  to  yield  reliable 
and  stable  values. 


TEST  SETUP 

The  test  setup  is  shown  in  Figures  3  and  4.  Figure  5  is  a  photo¬ 
graph  of  an  example  specimen.  Figure  3  shows  the  beams  were  tested  with 
the  notch  on  the  top  surface.  The  beam  weight  was  supported  by  four 
springs  aligned  with  reaction  rollers.  The  rollers  were  located  on 
bearing  pads  to  minimize  energy  dissipation  at  the  bearing  points. 

The  load  was  applied  through  a  closed-loop,  servo-controlled, 
20  kip,  MTS  testing  machine.  The  tests  were  displacement  controlled 
with  a  cross-head  displacement  rate  of  approximately  5.10  m/sec 
(2.10  ^  in. /sec) . 
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Two  clip  gages  with  a  sensitivity  of  0.0025  mm  (0.0001  in.)  were 
installed  across  from  the  load  point.  The  clip  gages  bore  against  two 
aluminum  beams  which  span  across  the  reaction  points  (Figure  5).  Tight 
fitting  slots  and  holes  machined  in  the  frame  allowed  for  rotation  and 
horizontal  displacement  without  vertical  movement  and  negligible  fric¬ 


tion. 


In  order  for  the  forces  on  the  beam  to  be  statically  determinate, 
the  two  reaction  rollers  bore  on  a  cylindrical  surface  whose  axis  was 
perpendicular  to  the  rollers'  axis.  Thus,  obtaining  a  single  point  of 
contact  equivalent  to  spherical  bearing. 

An  important  advantage  of  the  setup  lies  in  obtaining  the  long  tail 
of  the  LLPD  plot  (Figure  6),  representing  not  only  the  beam  weight  but 
also  other  similar  effects  (such  as  clip  gages  weight  and  reaction)  as 
well  as  possible  constant  friction.  If  some  variable  effects  are  pre¬ 
sent  during  the  test,  these  can  be  evaluated  by  observing  the  linearity 
of  the  plot  after  the  load  carrying  capacity  of  the  beam  has  been  spent. 


TEST  SERIES 


Three  series  of  tests  were  conducted.  RILEM  guidelines  were  fol¬ 
lowed  in  all  tests  for  maximum  aggregate  size,  conditions  of  storage, 
support  and  loading  arrangement,  accuracy  of  measurement,  and  rate  of 
loading. 

The  concrete  mixes  with  their  mechanical  properties  are  given  in 
Table  1.  The  maximum  aggregate  size  was  10  mm  (3/8  in.)  in  all  cases. 
The  initial  modulus  of  elasticity  in  compression,  E,  is  also  tabulated. 


Series  I 


Twelve  baseline  beams  (similar  to  RILEM's  smallest  specimen)  were 
cast  with  dimensions  of  102  by  102  by  788  by  838  mm  (4  by  4  by  31  by 
33  in.)  (width  by  depth  by  span  by  length)  (Figures  1  and  3).  The 


notch-to-depth  ratio,  aQ/W,  was  0.5.  During  the  tests,  the  beams  were 
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Table  1.  Material  Properties 


Series 

Cement 

(kg/m3) 

Water 

(kg/m3) 

@1 

10  mm  Gravel 
(kg/m3) 

Sand 

(kg/m3) 

f  ' 
c 

(MPa) 

ft 

(MPa) 

E 

(GPa) 

I 

279 

167 

0.60 

1062 

907 

29.0 

3. 1 

21.7 

II 

613 

245 

0.40 

1034 

443 

58.9 

4.2 

24.5 

III 

400 

220 

0.55 

1044 

_ 

540 

33. 1 

3.5 

19.7 

The  compressive  strength,  f  ,  was  obtained  at  28,  35,  and  30  days, 
respectively  for  each  series.  In  every  case,  three  152-  by  305-mm 
(6-  by  12-in.)  cylinders  were  tested. 

The  tensile  strength,  f.,  was  obtained  at  28,  32,  and  29  days, 
respectively,  using  the  same  type  of  cylinders.  Six  splitting 
tensile  tests  were  performed  for  Series  I,  then  only  three  each 
for  the  other  series  due  to  the  uniformity  of  the  values. 

The  modulus  of  elasticity  was  calculated  by  measuring  the  cylinders 
strain  at  the  beginning  of  the  compression  tests. 


unloaded  past  peak  load  and  both  and  G^*  obtained.  The  beams  were 

tested  in  four  groups  of  three  after  curing  for  27,  28,  29,  and  32  days 
in  a  fog  room. 

Series  II 

Twelve  additional  specimens  with  the  same  overall  dimensions  as 

Series  I  were  prepared.  They  were  divided  into  three  groups  of  four 

beams  with  the  same  a  /W.  The  Initial  notch-to-depth  ratio  was  0.3, 

o 

0.5,  and  0.7  for  groups  1,  2  and  3,  respectively.  In  order  to  isolate 
the  effect  of  notch  depth  from  the  effect  of  curing  time  one  specimen  of 
the  four  from  each  group  was  tested  on  the  same  day  after  curing  for  28, 
29,  32,  and  33  days. 


The  third  series  consisted  of  four  additional  baseline  specimens 
(same  dimensions  as  series  I  with  a  / W  =  0.5)  plus  four  specimens  with 
depths  of  216  mm  (8.5  in.)  but  with  the  same  width,  span,  length,  and 
aQ/W.  While  not  geometrically  similar  to  the  standard  RILEM  specimen, 
the  'crease  in  ligament  area  (unbroken  area  at  the  notch)  demonstrated 
the  size  effect  on  G^. 

One  baseline  specimen  and  one  deeper  beam  were  tested  each  day 
after  curing  for  28,  29,  30,  and  32  days. 

ANALYTICAL  REPRESENTATION  WITH  FINITE  ELEMENTS 

The  finite  element  program,  ADINA  (Ref  10)  features  user-supplied 
loading  which  can  be  expressed  as  an  arbitrary  function  of  nodal  dis¬ 
placements.  This  allowed  the  implementation  of  an  FCM  approach.  The 
element  mesh  was  derived  from  Reference  11  and  is  shown  in  Figure  7. 
Due  to  symmetry,  only  half  of  a  beam  needed  to  be  discretized.  The 
notch-to-depth  ratio,  aQ/W,  was  chosen  as  0.5.  The  beam  was  loaded  in 
displacement  control  (both  in  the  program  and  in  the  actual  tests)  and 
whenever  the  tensile  strength,  f^,  was  reached  at  a  node,  the  node  was 
released  and  the  midspan  displacement  was  step  increased  until  another 
fj.  was  exceeded  at  another  node.  This  iterative  process  was  continued 
until  the  crack  progressed  across  the  beam  cross  section. 

The  modulus  of  elasticity,  E,  used  in  the  finite  element  analyses 
was  measured  as  the  initiation  slope  of  the  stress  strain  diagram  of  a 
compression  test  on  a  152  -  by  305-mm  (  6-  by  12-in.)  cylinder.  Others, 
including  Peterson  (Ref  12),  have  used  the  dynamic  modulus  of  elas¬ 
ticity,  E^,  for  analysis.  The  difference  between  dynamic  and  static  is 
small  (E^  is  about  10  to  20  percent  higher  than  E)  and  the  loading  rates 
between  the  cylinder  tests  and  the  three-point  bend  (3PB)  tests  is  with¬ 
in  one  order  of  magnitude.  Small  differences  have  also  been  reported 
between  tests  carried  out  in  tension  compared  to  compression  (around  10 
percent  lower  according  to  Reference  13).  Hence,  E  was  chosen  as  the 
modulus  of  elasticity  of  the  equivalent  homogeneous  elastic  material. 


Four  stress  versus  crack  width  (o  versus  w)  relations  were  analyzed 
and  are  described  in  Figure  8.  The  straight  line  (SL)  and  bilinear  or 
concrete  (C)  relations  were  proposed  by  Peterson  (Ref  9),  the  exponen¬ 
tial  (E)  relation  is  proposed  by  the  authors  (from  curve  fitting  rela¬ 
tionships  derived  in  Ref  11  and  14),  and  the  power  function  (P)  by 
Reinhardt  (Ref  15).  In  all  cases,  the  relations  conform  to: 

Gf  =  f-a  dw 

as  recommended  by  Hillerborg  et  al.  (Ref  1). 

RESULTS 

Test  Series  I 

Values  of  G^*,  peak  load,  deflection  at  peak  load  (dp),  and  deflec¬ 
tion  when  the  load  carrying  capacity  of  the  beam  vanishes  (dQ),  are  re¬ 
ported  in  Table  2.  The  Increase  in  crack  length,  Aa  (measured  by  dye 
insertion)  is  also  Indicated.  In  some  cases  the  closed  loop  servo  con¬ 
trol  on  the  testing  machine  did  not  allow  the  unloading  to  take  place 
Immediately  after  the  maximum  load  (when  it  decreases  to  about  95  per¬ 
cent  of  its  peak  value),  so  experimental  values  were  not  reported. 

The  LLPD  plots  were  drawn  until  a  long  horizontal  trace  was  ob¬ 
tained  (Figure  6),  indicating  that  the  beam  was  not  carrying  any  load. 
The  noise  and  random  friction  then  measured  showed  an  effect  of  about 
±4  N  (±1  lb). 

Dye  insertion  highlighted  a  crack  front  which  appeared  fairly 
straight.  Typically  the  dye  was  inserted  from  the  top,  but  in  a  couple 
of  cases  excess  dye  was  allowed  to  run  along  the  sides  of  the  beam.  In 
those  cases  the  crack  front  appeared  curved  due  to  dye  absorption  on  the 
sides. 
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Test  Series  II 

Series  II  results  are  detailed  in  Table  3.  A  strong  dependency  of 

on  ao/W  is  apparent  with  variations  of  about  25  percent  (from  Its 

value  at  a  /W  =  0.5).  decreases  as  a  /W  increases.  It  should  also 
o  r  o 

be  noted  that  G^*  follows  G^  closely. 

Test  Series  III 

Series  III  results  are  reported  in  Table  4.  An  increase  in  G^  of 
about  24  percent  is  observed  when  the  depth  is  increased  from  102  to  216 
mm  (4  to  8.5  in.). 

Finite  Element  Analysis 

Figure  9  shows  the  average  LLPD  plot  from  the  12  tests  of  Series  I 
(solid  trace)  and  the  approximations  from  finite  element  analyses  car¬ 
ried  out  using  the  four  different  a  versus  w  relations. 


DISCUSSION 

Effects  on 

The  variation  of  G^  with  notch  depth  and  beam  depth  indicates  that 
the  fracture  energy  is  not  a  material  property  and  is  dependent  on  the 
specimen  configuration.  When  modelling  concrete  as  an  equivalent  homo¬ 
geneous  linear  elastic  material,  all  the  energy  supplied  is  assumed  to 
be  converted  into  surface  energy  by  propagation  of  a  single  macrocrack. 
The  microcracking  in  the  process  zone  also  dissipates  energy  that  the 
model  cannot  take  into  account.  The  process  zone  depends  upon  the 
stress  field  which  is  dependent  on  the  geometry  (Ref  16).  The  process 
zone  depends  upon  the  aggregate-to-specimen  size  ratio  which  Is  a  size 
parameter  (Ref  4).  This  dependency  raises  a  question  about  the  reli¬ 
ability  of  using  a  specimen  in  three-point  bending  instead  of  in  direct 
tension  to  determine  the  stress  versus  crack  width  relation. 
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Table  4.  Results  -  Series  III 


N 

(mm) 

Specimen 

Gf 

(N/m) 

V 

(N/m) 

Peak  load 
(N) 

Unloading 
(%  of  peak 

at 

1  oad) 

d 

P 

(mm) 

d 

0 

(mm) 

Aa 

(mm) 

1 

76.6 

76.8 

909 

90 

0.14 

2.9 

7.9 

3 

81.9 

- 

958 

- 

0. 16 

2.9 

- 

102 

5 

71.6 

70.6 

908 

89 

0. 15 

2.9 

7.9 

7 

59.2 

58.2 

793 

92 

0.12 

3.0 

3.2 

Mean 

72.3 

68.5 

892 

90 

0. 14 

2.9 

6.3 

2 

84.7 

- 

3795 

- 

0.09 

1.9 

- 

5 

89.5 

- 

3900 

- 

0.  10 

2.2 

- 

216 

8 

90.7 

- 

3875 

- 

0.11 

1.9 

- 

11 

92.2 

85.7 

3580 

89 

0.10 

2.0 

7.2 

— 

Mean 

89.3 

85.7 

3788 

89 

0. 10 

2.0 

7.2 

Sensitivity  analyses  (Ref  17)  and  further  analyses  by  the  authors 
have  shown  that  the  model  is  typically  not  as  sensitive  to  variations  in 
as  it  is  to  variations  of  the  other  parameters.  Most  of  the  energy 
does  appear  to  dissipate  through  crack  surface  formation  ensuring  the 
validity  of  the  model. 

versus  G^* 

Gf*  was  defined  in  an  attempt  to  quantify  and  eliminate  the  energy 
dissipated  outside  of  the  crack  zone.  The  strong  variation  of  G^  with 
aQ/W  seems  to  indicate  a  dependency  on  the  stress  field,  i.e.,  an  energy 
dissipation  through  microcracking  in  the  process  zone  which  is  greater 
than  indicated  by  G^*.  For  standard  specimen  sizes,  the  difference  be¬ 
tween  Gf*  and  G^  is  about  5  percent  while  the  fracture  energy  shows 
variations  in  the  order  of  25  percent.  Furthermore,  both  values  vary  in 
the  same  manner  and  show  size  and  geometry  dependency.  The  proposed 


3$$ 


approach  does  not  seem  to  discard  all  the  energy  spent  in  the  process 
zone.  Further  research  is  needed  to  refine  the  model.  Given  the  ad¬ 
ditional  steps  involved  in  its  evaluation,  G^*  does  not  appear  to  pre¬ 

sent  any  practical  advantage  over  G^. 

Crack  Front 

The  crack  front  appearance  seems  to  be  linked  to  the  dye  impregna¬ 
tion  procedure.  When  the  dye  was  inserted  only  from  the  top,  the  crack 
front  appeared  straight.  In  a  few  cases  the  dye  was  inserted  from  the 

sides  and  the  top,  and  then  the  front  appeared  curved  (Ref  16).  The 

latter  crack  front  appearance  seemed  to  indicate  an  anisotropy  of  the 
process  zone  where  the  microcracks  would  tend  to  merge  laterally  before 
joining  the  macrocrack.  In  the  determination  of  G^*  a  straight  front 
was  assumed. 

Analytical  Model 

The  good  agreement  between  the  ascending  part  of  the  LLPD  curve  of 
the  average  experimental  plot  and  the  model  results  supports  the  accept¬ 
ability  of  using  E  from  standard  cylinder  tests.  This  is  an  advantage 
since  E  is  more  readily  available  than  E^. 

Among  the  different  o  versus  w  relations,  the  SL  is  the  coarsest, 
11  owed  by  the  E  and  the  P  models.  The  bilinear  relation,  C,  gives 
only  a  slight  overestimation  of  the  peak  load  and  reasonably  approaches 
the  descending  branch.  Keeping  f^  and  G^  constant,  the  bilinear  rela¬ 
tionship  could  be  improved  to  yield  a  better  match,  as  done  by  Roelfstra 
and  Wittmann  (Ref  18)  and  Carpinteri  et  al.  (Ref  19). 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  fracture  for  notched  beams  in  three-point  bending  was  experi¬ 
mentally  determined.  Three  series  of  tests  were  conducted  on  a  total  of 
32  beams  and  a  finite  element  analysis  carried  out  using  a  nonlinear 
model.  The  following  conclusions  were  derived: 


•  The  experimental  setup  led  to  more  reliable  and  consistent 
values  of  concrete  fracture  energy  and  measured  crack,  length. 

•  G^*  (calculated  by  current  crack  size  and  ligament  area)  did  not 
present  any  substantial  advantage  over  (determined  by  RILEM 
guidel i nes) . 

•  Finite  element  analysis  yielded  accurate  representations  of  the 
fracture  behavior  when  coupled  with  a  bilinear  stress  versus 
crack  width  relation. 

•  Gf  1 s  dependent  upon  size  and  geometry.  The  fracture  energy 
increases  with  increasing  size. 

The  inconsistency  that  arises  from  applying  the  FCM  approach  seems 
to  be  linked  to  its  inability  to  represent  the  nonlinearity  of  the  ma¬ 
terial  and  the  energy  dissipation  outside  the  fracture  zone.  The  model 
assumes  all  the  energy  to  transform  into  surface  energy  through  forma¬ 
tion  of  a  single  major  crack  and  only  introduces  nonlinearity  in  the 
stress-crack  width  relation. 

The  fracture  model  for  plain  concrete  should  be  extended  to  three- 
dimensional  applications  and  to  bar-reinforced  concrete.  To  add  rein¬ 
forcement,  the  mechanism  for  transferring  shear  from  concrete  to  steel 
will  need  to  be  modeled,  but  the  technological  risk  is  low. 

ACKNOWLEDGMENT 

The  authors  wish  to  thank  Mr.  T.J.  Holland  for  the  assistance  pro¬ 
vided  in  the  implementation  of  the  numerical  calculations. 


REFERENCES 

1.  A.  Hillerborg,  M.  Modeer,  and  P.E.  Peterson.  "Analysis  of  crack 
formation  and  crack  growth  in  concrete  by  means  of  fracture  mechanics 
and  finite  elements,"  Cement  and  Concrete  Research,  vol  6,  1976,  pp 
773-782. 

2.  Y.R.  Rashid.  "Analysis  of  prestressed  concrete  pressure  vessels," 
Nuclear  Engineering  and  Design,  vol  7,  no.  4,  1968,  pp  334-344. 

3.  Z.P.  Bazant  and  L.  Cedolin.  "Blunt  crack  band  propagation  in  finite 
element  analysis,"  Journal  of  the  Engineering  Mechanics  Division,  Amer¬ 
ican  Society  of  Civil  Engineers,  vol  105,  no.  EM2,  1979,  pp  297-315. 

4.  Z.P.  Bazant  and  B.H.  Oh.  "Crack  band  theory  for  fracture  of  con¬ 
crete,"  Materials  and  Structures,  vol  16,  no.  93,  May-Jun  1983,  pp 
155-177. 

5.  Z.P.  Bazant,  J.K.  Kim,  and  P.  Pfeiffer.  "Continuum  model  for  pro¬ 
gressive  cracking  and  identification  of  nonlinear  fracture  parameters," 
Applications  of  Fracture  Mechanics  to  Cementitious  Composites 
(S.P.  Shah,  editor),  NATO-ARW,  Northwestern  University,  Dordrecht,  The 
Netherlands,  Sep  1984,  pp  197-246. 

6.  Y.  Jenq  and  S.P.  Shah.  "Two  parameter  fracture  model  for  concrete," 
Journal  of  Engineering  Mechanics,  American  Society  of  Civil  Engineers, 
vol  111,  no.  10,  Oct  1985,  pp  1227-1261. 

7.  RILEM  Technical  Committee  50-FMC.  "Draft  recommendation:  Deter¬ 
mination  of  the  fracture  energy  of  mortar  and  concrete  by  means  of 
three-point  bend  tests  on  notched  beams,"  Materials  and  Structures, 
no.  106,  Jul-Aug  1985,  pp  285-290. 


i 


aJvV- 

r-y.y. 


i 


v  »  % 

V 

VA> 

S,*V 
V  V, 

Vv 

W 


V 

>.V\j 
•  ’>  v 
Vv 


8.  S.  Mindess.  "Rate  of  loading  effects  on  the  fracture  of  cementi¬ 
tious  materials,"  Applications  of  Fracture  Mechanics  to  Cementitious 
Composites  (S.P.  Shah,  editor),  NATO-ARW,  Northwestern  University, 
Martinus  Ni.ihoff  publishers,  Dordrecht,  The  Netherlands,  Sep  1984,  pp 
617-636. 

9.  Kansas  State  University.  Research  Report  181:  Evaluation  of  pro¬ 
posed  methods  to  determine  fracture  parameters  for  concrete  in  bending, 
by  S.E.  Swartz  and  S.T.  Yap.  Manhattan,  KS,  Jun  1986. 

10.  ADINA  Engineering  Incorporated.  AOINA:  A  finite  element  program 
for  automatic  dynamic  incremental  nonlinear  analysis.  Watertown,  MA, 
Dec  1985. 

11.  Lund  Institute  of  Technology.  Report  TVBM- 100 1 :  A  fracture  me¬ 
chanics  approach  to  failure  analysis  of  concrete  materials,  by  M.A. 
Modeer.  Division  of  Building  Materials,  S-221  00,  Lund,  Sweden,  1979. 

12.  Lund  Institute  of  Technology.  Report  TVBM-1006:  Crack  growth  and 
development  of  fracture  zones  In  plain  concrete  and  similar  materials, 
by  P.E.  Peterson.  Division  of  Building  Materials,  S-221  00,  Lund, 
Sweden,  1981. 

13.  V.S.  Gopalaratnam  and  S.P.  Shah.  "Softening  response  of  plain 
concrete  in  direct  tension,"  Journal  of  the  American  Concrete  Institute, 
vol  82,  no.  3,  May-June  1985,  pp  310-323. 

14.  P.  Nallathambi  and  B.L.  Karihaloo.  "Prediction  of  load  deflection 
behavior  of  plain  concrete  from  fracture  energy,"  Cement  and  Concrete 
Research,  vol  16,  1986,  pp  7-16. 

15.  H.W.  Reinhardt.  "Fracture  mechanics  of  an  elastic  softening  mate¬ 
rial  like  concrete,"  Delft  University  of  Technology,  The  Netherlands, 
HERON,  vol  29,  no.  2,  1984,  pp  1-42. 


16.  A.  Bascoul,  F.  Kharchi,  and  J.C.  Maso.  "Concerning  the  measurement 
of  the  fracture  energy  of  a  micro-concrete  according  to  the  crack  growth 
in  a  three  point  bending  test  on  notched  beams,"  Proceedings  of  the 
International  Conference  on  Fracture  of  Concrete  and  Rock,  Society  for 
Experimental  Mechanics,  in  Houston,  Texas,  June  1987,  pp  631-643. 

17.  A.  Hillerborg.  "The  theoretical  basis  of  a  method  to  determine  the 
fracture  energy  of  concrete,"  Materials  and  Structures,  no.  106,  Jul- 
Aug  1985,  pp  291-296. 

18.  P.E.  Roelfstra  and  F.H.  Wittmann.  "Numerical  method  to  link  strain 
softening  with  failure  of  concrete,"  in  Proceedings  of  the  International 
Conference  on  Fracture  Mechanics  of  Concrete,  Lausanne,  Switzerland,  Oct 
1985,  pp  163-175. 

19.  A.  Carpinteri,  G.  Colombo,  G.  Ferrara,  and  G.  Giuseppettl.  "Numer¬ 
ical  simulation  of  concrete  fracture  through  a  bilinear  softening 
stress-crack  opening  displacement  law,"  in  Proceedings  of  the  Inter¬ 
national  Conference  on  Fracture  of  Concrete  and  Rock,  Society  for 
Experimental  Mechanics,  Houston,  Texas,  June  1987  pp,  178-191. 


DISTRIBUTION  LIST 


AF  (’550  (  I  S  Hi  l  i  Patrick  AFB.  FI. 

AFESC  1ST  (Library),  Tyndall  AFB.  FI 

ARMY  AMCS.M  WS.  Alexandria.  VA:  ERADCOM  lech  Supp  Dir  (DELSD-1).  Ft  Monmouth.  N.l 
ARMY  BELVOIR  RieD  CEN  STRBE-AALO.  Ft  Bchoir.  VA:  STRBE-BLORE.  Ft  Belvoir.  VA 
ARMY  CERL  Library.  Champaign.  II. 

ARMY  CRREL  Library.  Hanover.  NH 

ARMY  ENGR  DIST  Phila  I  T.  Philadelphia.  PA 

ARMY  MISSILE  R&D  CMD  Ch.  Docs.  Set  lido  Or.  Arsenal.  AL 

ADMINSL'PU  PWO.  Bahrain 

BUREAU  OF  RECLAMATION  D- 1512  (CIS  DePuy).  Denver.  CO;  Smoak.  Denver.  CO 

CBC  Library.  Davisville.  RI 

CBL  411.  OIC.  Norfolk.  VA 

COGARD  RiVDC  Library.  Groton.  CT 

COMDT  COGARD  Library.  Washington.  DC 

NAVRESCEN  PE-PLS.  Tampa.  FL 

D1A  DB-6E1.  Washington.  DC;  DB-6E2.  Washington.  DC:  VP-TPO.  Washington.  DC 

D1RSSP  Tech  Lib.  Washington.  DC 

DNA  STTITL.  Washington.  DC 

DTIC  Alexandria.  VA 

DTRCEN  Code  4111.  Bethesda.  MD 

FCTC  LANT.  PWO.  Virginia  Bch.  VA 

GIDEP  OIC.  Corona.  CA 

LIBRARY  OF  CONGRESS  Sci  &  Tech  Dtv.  Washington.  DC 
MARCORBASE  Code  4.01 .  Camp  Pendleton.  CA 

NAS  Code  165.  Keflavik.  Iceland:  Miramar.  Code  IX2IA.  San  Diego.  CA:  PWO  (Code  IN2)  Bermuda:  PWO. 
Key  West.  FL;  PWO.  Sigonella.  Italy 

NATL  BUREAU  OF  STANDARDS  Bldg  Mat  Div  (Malhey).  Gaithersburg.  MD 
NAVCAMS  SCE  (Code  N-7).  Naples.  Italy 
N AVCOASTSYSCEN  lech  Library.  Panama  City.  FI, 

NAVCOMMSTA  Code  401.  Nea  Makri.  Greece 
NAVEODTECHC  EN  Tech  Library  .  Indian  Head.  MD 

NAVF ACENGCOM  Code  05.  Alexandria.  VA:  Code  04A.  Alexandria.  VA:  Code  (MAID.  Alexandria.  VA: 
Code  04A5.  Alexandria.  VA;  Code  (I4A4E  (Bloom).  Alexandria.  VA;  Code  04B2  (M.  Yaehnis). 

Alexandria.  VA;  Code  0651  (Cyphers).  Alexandria.  VA;  Code  I002B.  Alexandria.  VA:  Code  1115. 
Alexandria.  VA 

NAM  ACENGCOM  ■  CUES  DIV  FPO-IPL.  Washington.  DC 
NAVFAC ENGCOM  LANT  DIV  Library.  Norfolk.  VA 
NAVFAC T.NGCOM  -  NORTH  DIV  Code  (UAL.  Philadelphia.  PA 
NAVFACI  NCiCOM  ■  PAC  DIV.  Library.  Pearl  Harbor.  Ili 
NAM  AC  l  NC.COM  -  SOUTH  DIV  Library  .  Charleston.  SC 
NAM. ACENGCOM  -  WEST  DIV  C  ode  04A2.2  (Lib).  San  Bruno.  C  A 
N  \\ FAC  I  NC.COM  CONTRACTS  OIC  C  .  Rota.  Spam 
ARMS  C  ORPS  OF  I  NC  IRS  Library.  Seattle.  WA 
N  A  1 1  BURI  AL  OF  STANDARDS  R  Chung.  Gaithersburg.  MD 
NAVC  T I A  PC  i  Rl  Code  Ml.  Williamsburg.  VA 
I  Ieoh  Lib.  Pensacola.  FT. 

N AVFACENC iCC )M  Corle  00.  Alexaiulria.  VA:  Code  051  (L.ssoglou).  Alexandria.  VA;  Code  04AI. 
Alexandria.  VA:  C  ode  ICA  (Herrmann).  Alexandria.  VA:  C  ode  07M  (Gross).  Alexandria.  VA:  Code 
IIUMI24  (lib).  Alexaiulria.  \  A 

NAVFAC T  N'CICC )M  -  I  AN!  DIV  Br  Ole.  Dir.  Naples.  Italy 

NAVFAC  ENGCOM  CON  I  RAC  I  S  DROIC  C  .  I.emoore.  (  A:  ROIC  C  .  \  irginia  Beach.  VA 
NAVOCT  ANSYSC  I  N  DEI  .  lech  Lib.  Kailua.  Ill 
N  A\  SII1PRI  PFAC  I.ihran.  Guam 

NMSIllP'i  D  C  ode  202,4.  l  ong  Beach.  C  A;  Library.  Portsmouth.  NH 
NAVSTA  CO.  Rooserelt  Roads.  PR;  Engrg  Dir.  Rota.  Spam 
N  A\  SUPPAC  T  PWO.  Naples.  Italy 
NAVSUPPC)  Sec  Ollr.  La  Maddalena.  Italy 

N.AVSWO  Code  L2II  (Miller).  D.ihlgrcn.  VA:  Code  W  42  (R  Pon/etlol.  Dahlgren,  \  A.  PWO.  Dahlgien.  \  A 

NRI.  Code  5 NIKI.  Washington.  DC 

OC  NR  DM.  Code  4X1.  Bat  St  loins.  MS 

PAC  MISRANFAC  III  Area.  PWO.  Kekaha.  Ill 

PERRY  OC  EAN  ENGRG  R  I’ellcn.  Riviera  Beach.  I  I 

PM  I C  C  ode  HUN.  Poult  Mill’ll.  CA 


PWC  ACE  Office.  Norfolk.  VA:  Code  101  (library).  Oakland.  (  A:  Code  |M-(  .  San  Diego.  C  A:  (  ode  On. 
Circal  Lakes.  II.:  Library  (Cotie  134).  Pearl  Harbor.  Ill;  l.ibiarv.  Ouaiii.  Mariana  Mantis:  I  ibian.  Norlolk. 
VA:  Libran.  Pensacola.  FL:  l.ibrarv.  Yokosuka.  Japan:  lecli  librarv.  Subic  Ua\.  RP 
I  S  DL PI  OF  INTERIOR  Natl  Park  Sve.  RMR  PC.  Pemer.  CO 
CSS  CSS  FIT  .TON.  Code  W-3 

CAI  IFORNIA  Nav  &  Ocean  Dev  (Armstrong).  Sacramento.  CA 
Cl  ARKSON  COLL  OF  TECH  CL  Dept  I  Batson).  Potsdam.  NY 
COLORADO  SC  HOOL  OF  MINES  Dept  of  Engrg  (Cluing).  Golden.  CO 
COLORADO  STATE  CNIVERSITY  CE  Dept  <W  Charlie).  Fort  Collins.  Ml) 

CORNELL  UNIVERSITY'  Civil  &  Environ  Engrg  (Dr.  Kulhavvv).  Ithaca.  NS  :  l.ibrarv.  Ithaca.  NY 
DI  KE  CNIVERSITY  CE  Dept  (Muga).  Durham.  NC 

FLORIDA  ATLANTIC  UNIVERSITY’  Ocean  Engrg  Dept  (Hart).  Boca  Raton.  FT  :  Ocean  Engrg  Dept 
(McAllister).  Boca  Raton.  FI.:  Ocean  Engrg  Dept  ( Su ) .  Boca  Raton.  FL 
FLORIDA  INST  OF  TECH  CE  Dept  (Kalajian).  Melbourne.  IT.;  .1  Schwalbe.  Melbourne.  FL 
JOHNS  HOPKINS  UNIV  C'E  Dept  (Jones).  Baltimore.  MD;  Cites  Bav  Rsch  Inst.  Rsch  Lib.  Shatlv  Side.  MD 
LAWRENCE  LIVERMORE  NATL  LAB  FJ  Tokarz.  Livermore.  CA:  L-(o4.  Plant  Engrg  Lib.  Livermore.  CA 
LEHIGH  UNIVERSITY  Linderman  Library.  Bethlehem.  PA 
MIT  Engrg  Lib.  Cambridge.  MA 

NAI  L  ACADEMY  OF  SC  IENCES  NRC  .  Naval  Studies  Bd.  Washington.  DC 
OKLAHOMA  STATE  UNIV  CE  Seol  (Lloyd).  Stillwater.  OK 
PENNSYLVANIA  STATE  CNIVERSITY  Applied  Rsch  Lab.  State  C  ollege.  PA 
Pl’RDCE  UNIVERSITY  Engrg  Lib.  W  Lafayette.  IN 
STATE  UNIVERSITY  OF  NEW  YORK  CE  Dept.  Buffalo.  NY 
TEXAS  Act  I  UNIVERSITY  Civil  &  Mech  Engr  Dept.  Kingsville.  IX 
TEXAS  \cV M  CNI\ ERSITY’  Ocean  Engr  Proj.  College  Station.  I  X 
l  Nl\'  OF  TENNESSEE  CE  Dept  (Kane).  Knoxville.  IN 
UNIVERSITY  OF  C  ALIFORNIA  Engrg  Lib..  Berkeley.  CA 

UNIVERSITY’  OF  DELAWARE  C'E  Dept.  Ocean  Engrg  (Dalrvmplc).  Newark.  DL:  Engrg  Col  (Dexter). 
Lewes.  DE 

CNIVERSITY  OF  HAWAII  Manoa.  l.ibrarv.  Honolulu.  III 
l  NIVERSITY  OF  II  I.INOIS  l.ibrarv.  Crbana.  11  .  Met/  Rel  Rm.  f  rbana.  II. 

I  NIVERSITY  OF  RHODE  ISLAND  CE  Dept  (Kovaes).  Kingston.  Rl 

l  ’  N 1 VE  RSIT  Y  OF  I’E.XAS  CL  Dept  (R  Olson).  Austin.  IX.  Cl  Dept  (Thompson).  Austin.  IX 
l  NIX  I  RSI  1 Y  Of  C  ALIFORNIA  C  l  Dept  (Fcitves).  Berkeley .  CA 
l  NIYERS1 1’Y  OF  II  X AS  ECJ  5.41 1:  (Friedrich).  Austin.  IX 
l  N | VI  RSI  I  Y  Of  W  ASHING  ION  CL  Dept  (N  Hawkins).  Seattle.  W  A 
WISII’RN  ARCHEOLOGICAL  (  I  N  l  ibrary.  Tucson.  AZ. 

\MI  R I (  AN  STS  I  NGRG  CORP  B  Williamson.  Virginia  Beach.  VA 

\M1  RICAN  (  T )N(  Ri  l  l  INST  ITT  I  F  l.ibrarv.  Detroit.  Ml 

H  X  I  II  I  LI  New  I  Tie  Marine  Rsch  lab.  lab.  Duxbury.  MA 

tail  |  |)  I N  C  (Tics  Instrn  Div .  lech  Lib.  Glen  Burnie.  MD 

I  IND  \  HA!  1  I  IBRARY  Doc  Dept  Kansas  City.  MO 

MW  /I  ALAND  NY  Concrete  R'di  Assoc.  Library.  Porirtia 

XI  All  (  II  CORP  Petoni.  Miami.  FL 

I  IDI  W  \  1 1  R  CONSTR  CO  .1  Fowler.  Virginia  Beach.  VA 
I  RW  I N(  |  nar  I  ibr.iry.  Cleveland.  OH 

WISS.  JANNE’Y.  ITSINER.  X  ASSOC  DW  Pleiler.  Northbrook.  II 
W  ( )( >DW  ARD-CLYDI  CONSt  I  FANIS  R  Cross.  Walnut  (reek.  CA 
W|  SIT  OIL  WNI  Miami.  FT 


INSTRUCTIONS 


The  Naval  Civil  Engineering  Laboratory  has  revised  its  primary  distribution  lists.  The  bottom  of  the  label 
on  the  reverse  side  has  several  numbers  listed.  These  numbers  correspond  to  numbers  assigned  to  the  list  of 
Subject  Categories.  Numbers  on  the  label  corresponding  to  those  on  the  list  indicate  the  subject  category  and 
type  of  documents  you  are  presently  receiving.  If  you  are  satisfied,  throw  this  card  away  (or  file  it  for  later 
reference). 

If  you  want  to  change  what  you  are  presently  receiving. 

•  Delete  —  mark  off  number  on  bottom  of  label. 

•  Add  -  circle  number  on  list. 

•  Remove  my  name  from  all  your  lists  —  check  box  on  list. 

•  Charge  my  address  -  line  out  incorrect  line  and  write  in  correction  (PLEASE  ATTACH  LABEL). 

•  Number  of  copies  should  be  entered  after  the  title  of  the  subject  categories  you  select. 

Fold  on  line  below  and  drop  in  the  mail. 

Note:  Number*  on  label  but  not  listed  on  questionnaire  are  for  NCEL  use  only,  please  ignore  them. 


Fold  on  line  and  staple. 


DEPARTMENT  OP  THE  NAVY 


NAVAL  civil  engineering  laboratory 
PORT  HUENEME.  CALIFORNIA  93043-6003 


OFFICIAL  BUSINESS 
PENALTY  FOR  PRIVATE  USE.  *SOO 
NCEL-2700/4  (REV.  10-87) 
Oea0-U^L70~0044 


Commanding  Officer 
Code  L08B 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  California  93043-5003 


BUSINESS  REPLY  CARD 

FIRST  CLASS  PERMIT  NO  69 _ 

POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


DISTRIBUTION  QUESTIONNAIRE 

The  Naval  Civil  Engineering  Laboratory  is  revising  its  primary  distribution  lists. 


SUBJECT  CATEGORIES 

1  SHORE  FACILITIES 

2  Construction  methods  end  materials  (including  corrosion 

control,  coatings) 

3  Waterfront  structures  (maintenance/deterioration  control) 

4  Utilities  (including  power  conditioning) 

5  Explosives  safety 

6  Aviation  Engineering  Test  Facilities 

7  Fire  prevention  and  control 

8  Antenna  technology 

9  Structural  analysis  and  design  (including  numerical  and 

-  computer  techniques) 

10  Protective  construction  (including  hardened  shelters, 

shock  and  vibration  studies) 

1 1  Soil/rock  mechanics 

14  Airfields  and  pavements 

15  ADVANCED  BASE  AND  AMPHIBIOUS  FACILITIES 

16  Base  facilities  (including  shelters,  power  generation,  water  supplies) 

1 7  Expedient  roads/airfields/bridges 

18  Amphibious  operations  (including  breakwaters,  wave  forces) 

19  Over-the-Beach  operations  (including  containerization, 

materiel  transfer,  lighterage  and  cranes) 

20  POL  storage,  transfer  and  distribution 


28  ENERGY/POWER  GENERATION 

29  Thermal  conservation  (thermal  engineering  of  buildings,  HVAC 

systems,  energy  loss  measurement,  power  generation) 

30  Controls  and  electrical  conservation  (electrical  systems. 

energy  monitoring  and  control  systems) 

31  Fuel  flexibility  (liquid  fuels,  coal  utilization,  energy 

from  solid  waste) 

32  Alternate  energy  source  (geothermal  power,  photovoltaic 

power  systems,  solar  systems,  wind  systems,  energy  storage 
systems) 

33  Site  data  and  systems  integration  (energy  resource  data,  energy 

consumption  data,  integrating  energy  systems) 

34  ENVIRONMENTAL  PROTECTION 

35  Solid  waste  management 

36  Hazardous/toxic  materials  management 

37  Wastewater  management  and  sanitary  engineering 

38  Oil  pollution  removal  and  recovery 

39  Air  pollution 

44  OCEAN  ENGINEERING 

45  Seafloor  soils  and  foundations 

46  Seafloor  construction  systems  and  operations  (including 

diver  and  manipulator  tools) 

47  Undersea  structures  and  materials 

48  Anchors  and  moorings 

49  Undersea  power  systems,  electromechanical  cables, 

and  connectors 

50  Pressure  vessel  facilities 

51  Physical  environment  (including  site  surveying) 

52  Ocean. based  concrete  structures 
54  Undersea  cable  dynamics 


TYPES  OF  DOCUMENTS 

85  Techdau  Sheets  86  Technical  Reports  and  Technical  Notes 

83  Table  of  Contents  &  Index  to  TDS 


82  NCEL  Guides  &  Abstracts 
91  Physical  Security 


rj  None  — 

remove  my  name 


